Phonon deficit effect and solid state refrigerators based on superconducting tunnel 

junctions. 
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Since the first demonstration of electron cooling in superconductor-insulator-normal metal (SLN) 
tunnel junctions, there is growing interest in tunneling effects for the purpose to develop on-chip 
refrigerators which can generate cooling in the order of WOmK in an environment of 0.3-0.5K. Thin 
film devices have the advantage of being extremely compact, operate in a continuous mode, dissipate 
little power, and can easily be integrated in cryogenic detectors. Motivated by such possibilities, we 
investigate the phonon deficit effect in thin film SIS (superconductor-insulator-superconductor) 
and SIN tunnel junctions. Under certain circumstances, the phonon absorption spectra of such 
tunnel junctions have spectral windows of phonon absorption/emission. We propose to use phonon 
filters to select the phonon absorbtion windows and thus to enhance the cooling effect. Membranes 
attached to such tunnel junctions can be cooled in this way more effectively. We discuss a particular 
superlattice design of corresponding phonon filters. 



I. INTRODUCTION 



Solid state micro-refrigerators (SSMR) have many 
advantages over conventional cooling devices: they have 
a long life-time, they are compact and simple to ap- 
ply. These SSMRs utilize motion of electrons rather 
than atoms and molecules. The flow of electric cur- 
rent through a semiconductor-semiconductor contact is 
accompanied by a transfer of heat (the Peltier effect) 
from one semiconductor to the other. Micro-refrigerators 
based on the Peltier effect are working at tempera- 
tures above 150K . Solid-state micro-refrigerators op- 
erating below 150K are a scientific and technical chal- 
lenge. A possible approach to this problem is to 
consider electron kinetics in superconductor-insulator- 
normal metal (SIN) or superconductor-insulator- 
superconductor (SIS) tunnel junctions (|, As we 
will demonstrate in this article, another even deeper un- 
derstanding of the same effects is possible from the con- 
cept of the phonon deficit effect (PDE). The phonon 
deficit effect was initially predicted fi lOfl to occur when 
superconductivity is enhanced in high frequency electro- 
magnetic fields (UHF) M. Later the PDE was recog- 
nized to exist in more general cases O, |l2| . It is related 
to the violation of detailed balance of kinetic processes 
in superconductors in the presence of an external field. 
If a superconducting film is immersed into a heatbath 
and the electron system goes over to an excited nonequi- 
librium state due to an external supply of energy, then 
under appropriate conditions the superconducting film 
absorbs phonons from the heatbath rather than emitting 
phonons. 

The situation is essentially similar in the case of thin 
tunnel junctions of superconductors. Phonon emission 
of a symmetric SIS tunnel junction has been consid- 
ered in detail in Ref. ]IT[ and an asymmetric junction 
has been considered in Rcfs. |^|, [jj], Q. In the sub- 
threshold tunneling regime (V < 2A/e, where A denotes 
the gap) the resulting phonon flux is shown in Fig.|lJ. 
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FIG. 1: Spectral dependence of the phonon flux (in units of 
energy per unit time) radiated from a thin superconducting 
film (with smaller gap A s ) in an SIS' tunnel junction (taken 
from Ref. jjj)) 



One observes the phonon emission in the frequency win- 
dows < ojq < 2A/H and aj q > 3 A/ft. In the same time 
there is the phonon absorption in the frequency window 
2 A/h < bj q < 3 A/h. The analogy between this case and 
the superconductor in the UHF field [see Ref. |Q , Fig- 
ure 6.2] strikingly demonstrates the universality of the 
mechanism of the phonon deficit effect. 

Experiments on refrigeration are usually performed on 
the SIN tunnel junctions. So called "cooling by extrac- 
tion" is recognized as the operational principle of these 
type of refrigerators. At temperatures much below the 
critical temperature T c of the superconducting electrode 
(S) there are few thermally excited quasiparticles above 
the superconducting gap A, so that tunneling from S to 
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the normal metal electrode (N) can be neglected. Mean- 
while, in the normal metal electrode (N) there are enough 
electron-like and hole-like thermal excitations near the 
Fermi energy Ep which can tunnel to S and effectively 
modify the electron distribution in N. The rate of elec- 
tronic heat exchange with the lattice is proportional to 
Tg — Tp hl where T e is electronic temperature and T p h is 
the temperature of phonons. The cooling process is a bal- 
ance between phonon heat flow, quasiparticle extraction 
and Joule heating. The value of the applied potential V 
on the tunnel junction determines if cooling or heating 
occurs depending on which process is dominant in the 
heat transfer balance. 

Here we have studied the behavior of SIN tunnel 
junction using the phonon deficit effect. The outline of 
this paper is as follows. In Section || we describe the 
phonon heatbath model and decoupling of electrons from 
phonons. In Section III we solve kinetic equations for the 



SIN tunnel junction. We calculate the net phonon flux 
to the SIN tunnel junction when electron and phonon 
systems are maintened at different temperatures. We 
show that the rule for heat exchange P oc T^—T^ h follows 
from the PDE in SIN tunnel junctions. In analyzing the 
PDE we propose the phonon filters. We present the re- 
sults of our numerical simulations, which can be utilised 
to enhance the refrigeration effect. The constraction of 
phonon filters in accordance with the calculated phonon 
absorbtion spectra is discussed in Section |vj Conclu- 
sions are given in Section |y|. 



II. PHONON HEATBATH MODEL AND 
DECOUPLING OF ELECTRONS FROM 
PHONONS 

Experimentally junction electrodes are thin films 
placed on a substrate. These thin films should be consid- 
ered as a system of two interacting subsystems: electrons 
and phonons. Noncquilibrium dynamics of phonons and 
electrons in superconductors is described by a coupled 
system of kinetic equations Jl^| , In the simplest case, 
when the lattice phonons play the role of a thermostat, 
the kinetic equations for electron (hole) excitations and 
phonons in a metal can be written as in the following 
form 111 



dn± 



dt 

dNl 



Qf(n±e)+I el - ph (n ±e ,Nl q ), (1) 
i=L(Ni q )+P h -° l (N^,n ±e ) . (2) 



Here n± e denotes the distribution functions of nonequi- 
librium electron-like (n e ), respectively hole-like (n— e ) 
excitations in the film. N^ represents the distri- 
bution function of internal phonons, I el ~ ph (n± e , N^ ) 
stands for the electron-phonon collision integral, and 
jph-ei (JVj* , 7i± e ) denotes the phonon electron collision 

integral. Q> (n± e ) is an external source of nonequilib- 
rium quasiparticles and / represents external fields to 



be specified later. L(N^ ) is an operator describing the 
interaction of phonons with an external heatbath. 

An important problem in a dynamical theory of met- 
als is the study of noncquilibrium stationary states, i.e. 
when dn± € /dt = 0, dN^ /dt = 0. For thin films in an 
external field the dynamical problem can be simplified 
by considering the phonon heatbath model j^, |l2| and 
also the " geometric-acoustical" approximation of phonon 
propagation. In this approximation the distribution func- 
tion of internal phonons is given by N^ — N® , where 
iVj is the distribution function of phonons at thermo- 
dynamic equilibrium. The distribution of electrons and 
holes corresponds to a nonequilibrium state. In a steady 
state, using the phonon heatbath model, the following 
system of equations holds 



Q f {n±e) =I el - ph (n ±e ,N^, 



(3) 
(4) 



These equations are uncoupled. Eq.(||) can be used to 
obtain stationary steady state solutions for the distri- 
bution functions of nonequilibrium excitations (electron 
and hole branches. Eq.([|) becomes an identity, which de- 
scribes the phonon exchange (phonon flux) between the 
thin film and its environment (131 . 



III. SIN TUNNEL JUNCTION - KINETIC 
EQUATIONS IN THE STEADY STATE 

We consider now a massive superconductor (injector 
S'), having a superconducting gap A', joined with a ther- 
mostat held at temperature T. A thin film having a 
thickness d ~ £o (superconductor or normal metal) is at- 
tached to the massive superconductor by means of a thin 
oxide layer. It has the same temperature (see Fig.Q]). In 
the dissipative steady state, noncquilibrium electron (n e ) 
and hole (yi-e) excitations in the thin film (S or N) are 
described by the following kinetic equation (see Eqs. (|l|) 
and (0) ) 



(5) 



Here u e = \e\9(e 2 - A 2 )/Ve 2 - A 2 denotes the BCS den- 
sity of states, e is the quasiparticle energy and Q(n± e ) = 
Q ( n ±e) describes the injector. In a particular geometry 
the external field / represents the distribution function 
of quasiparticles in the electrode denoted by S'. The 
explicit form of the electron-phonon collision integral 
r l ~P h (n t ,N° q ) and Q(n ±e ) are given in Refs. § p). 
The dissipative steady state solutions of Equation (||) 
give the electron and the hole distribution functions in 
the thin film electrode. The corresponding phonon fluxes 
emitted to (or absorbed from) its environment in a fre- 
quency interval duj q centered at frequency Lu q can be ob- 
tained from (0). 
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The number of phonons absorbed per time unit in the 
volume f2 at frequency uj q in the spectral interval dto q is 
given by 

<W U , = -L (JV2.) PoK)d^ = Po (Lo q )P h - el (a* ) du q . 

(6) 

Here po(u> q ) = ilui q /(2ir 2 u 3 ), u is the velocity of sound 

and X = dX/dt. The absorbed power per volume f2 at 
frequency uj q in the spectral interval dtu q is given by 



dP 

duJn 



(7) 



Introducing the operator I(N® ) via 



allows to express the absorbed power (Q) in dimensionless 
form 



dP _ n\uj D {k B T d f 
dx WireFU 3 h 3 

n\uj D {k B T c ,f 



I(N° X ) = P I(N° X ) 



P n = 



167refu 3 h 3 



(8) 



where we abbreviate x = ftto q /(kBT c i) and T c < denotes 
the critical temperature of the injector. 

In the following we use the units fcs=ft=e=lif not 
stated otherwise. Also, recall that the primed variables 
correspond to the injector material. 



A. Numerical results and analysis 

To simplify further the analysis, we need to analyze 
properties of equation (^J). The electron-phonon col- 
lision integral I el ~P h is proportional to the parameter 
7 oc T 3 /uijj (damping) which characterizes the strength 
of the relaxation processes. The parameter v describes 
the intensity of the particle and hole injection Q(n± e ). 
The typical values of v are small compared to 7. When 
1//7 tends to zero, the particle and hole distribution func- 
tions go over to thermal equilibrium distribution func- 
tions and both, I el ~ ph (n± e , N° q ) and I ph - el (N° q , n± e ) 
vanish. At finite values of u/j the particle and hole dis- 
tribution functions differ from the thermal equilibrium 
functions and thus phonon fluxes arise. As observed in 
Refs. H ||, |l8| both negative and positive phonon fluxes 
occur between the thin film and the thermal reservoir. 
The net phonon flux integrated over all frequencies can 
be either negative or positive. 



V < A'(T). When the potential difference V across 
the tunnel junction approaches the superconducting gap 
A'(T), the phonon absorption decreases. The absorbed 
power increases with the increasing values of v/j- The 
situation changes dramatically in the high temperature 
regime T/T c i > 0.5. When the value of the potential 
applied across the tunnel junction approaches the su- 
perconducting gap A', then the S'lN tunnel junction 
starts to emit phonons to the environment. This is 
shown in Fig.|Qb]. In Fig. Q we present the dynamics of 
phonon flux in the vicinity of the minima of Figs |^a,b] . 
As expected from Fig.[||a], the phonon absorption spec- 
trum is spread out over the whole frequency range of 
acoustic phonons at the temperature T/T c i = 0.5 and 
A' — V > 0.18. At temperature T/T' c = 0.8 there is 
both phonon emission and phonon absorption at high and 
low frequency regions, respectively. The phonon emission 
spectrum of a SIN tunnel junction is different from the 
S'lS case (see Fig. [[!)). The SIS tunnel junction has 
two emission and one absorption frequency regions. In 
57 AT tunnel junctions the low frequency emission region 
is absent. Also, as temperature decreases the phonon 
emission decreases faster than the absorption. 



B. PDE in the SIN tunnel junctions:" Equilibrium" 

limit 

An analytic expression for the power absorbed by the 
57 tunnel junction when electron and phonon subsys- 
tems are maintained at temperature T e and T p h, respec- 
tively, can be obtained by integrating expression (||) . The 
heatflow from the phonon subsystem to the electron sub- 
system maintained at different temperatures is given then 
by (for details see Appendix ) 



P = a 



(9) 



where a = APq T(5) C(5)> C( x ) denotes Riemann's Zeta 
function and T(n) is the Gamma function. A comparison 
of Equation (j^) with the Equation (19) of Ref. jjj shows 
that P cx T e 5 - Tp h is compatible with the PDE in 57iV 
tunnel junctions [see Refs. |l6|, p3|]. This result shows 
that the usual approach || |q, [15], |l6| takes into account 
only the net phonon flux. In this approach details like 
absorption at low frequency region and emission at high 
frequency region (see Fig. ) are inaccessible. However, 
the details of phonon absorption spectra are important 
and may open new possibilities for cooling by 57A tunnel 
junctions as discussed in the next Section. 



Figure M shows the power P(V)/P(0) absorbed by 
the normal metal thin film electrode. It is shown as a 
function of the potential applied across the S'lN tun- 
nel junction for different values of the parameter 1^/7. 
All curves have a minimum (maximal absorption) at 



IV. MICRO REFRIGERATOR BASED ON THE 
PHONON DEFICIT EFFECT 

Why is it important to know the spectral dependence 
of the phonon fluxes? How can the phonon-deficit-effect 



4 





FIG. 2: Absorbed power P(V), in units of Po, by the nor- 
mal metal electrode of the thin film as function of applied 
potential for different values of u/j. A (T) and T c ' are the 
superconducting gap and the critical temperature of the in- 
jector, respectively, uiu is the Debye frequency of the normal 
metal electrode of the thin film and V is the applied potential 
across the tunnel junction, (a) T/T c , = 0.5, (b) T/T c , = 0.8. 
The lines are a guide to the eye. 



influence the microrefrigerator design? The answers to 
these questions can be found in spectral- selective phonon 
filters. 

To clarify this point, we should calculate the energy, 
which the nonequilibrium phonons yield to the heatbath. 
This quantity is equal to the integral of the spectral func- 
tion plotted in Fig. H . This integral having a net negative 



FIG. 3: Spectral absorption (negative part) and emission 
(positive part) of phonons near the maximal phonon absorp- 
tion potential of the thin film. A'(T) and T c i are the super- 
conducting gap and the critical temperature of the injector, 
respectively, lud is the Debye frequency of the normal metal 
electrode of the trhin film and V is the applied potential across 
the tunnel junction, (a) T/T c , = 0.5, (b) T/T c , = 0.8. 

value is likely to be higher in the case of an asymmet- 
ric S'lS junction ^ Q and also in asymmetric S'lN 
junctions. Fig.[[lj displays the nonequilibrium phonon 
flux from the S'lS tunnel junction, As << A.gt. As 

follows from Figs.H], 0andJ01 the net flux integrated 
over energy is negative"]]?], ||, Q. To enhance the PDE 
effect and achieve better cooling, one should separate 
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nonequilibrium negative and positive fluxes across the 
interface between the superconductor and the attached 
heat reservoir. Let us consider placing a phonon spectral 
filter between superconductor S and reservoir 2 (Fig. 
Suppose we can construct a filter p3, 17 which is trans- 



Transmission spectra 
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Here we consider the simplest case of a phonon fil- 
ter (I?], [l8) for the case when acoustic phonons propa- 
gate through a superlattice perpendicular to interfaces. 
In Fig.|| the superlattice system is shown. The system 




FIG. 4: Layout of the microrefrigerator with the improved 
properties. Inserts show suggested spectral characteristics of 
superlattice filters. 

parent to phonons in a window of frequencies where the 
flux in Fig.jjlJ is negative. Such a filter would permit the 
flow of phonons (Su q « 2A) from the reservoir 2 to the 
superconductor S. At the same time the filter will ab- 
sorb phonons below and above that frequency window, 
including those coming from the heated barrier. Another 
filter with complementary (i.e. the band pass) transmis- 
sion properties may be placed between reservoir 1 and the 
electrode S' . It will prevent phonons coming from S' (i.e., 
the reservoir 1) from being absorbed in S. This will make 
the absorption from the reservoir 2 more efficient. Then 
the system shown in Fig.JlJ will work as a refrigerator, by 
coolingthe reservoir 2. This cooling effect cannot reduce 
the temperature much below the critical temperature of 
superconductor S' (Fig. [Q ) . However, there are differ- 
ent classes of superconductors with critical temperatures 
covering a wide range of temperatures from about 150 A" 
to very low temperatures. Then, cooling cascades can be 
organized in such a way to generate substantial cooling. 
We suggest that such designs of a PD_E-based microre- 
frigerator could be carefully considered for the purpose 
of practical use. In the next section we will discuss some 
designs. 



A. Design of phonon filters 



FIG. 5: Bragg multilayer between two media. 

consists of a sequence of different layers: material 1 with 
thickness d\ and impedance Z\ , material 2 with thickness 
d,2 and impedance Z2, etc. If one repeats this configu- 
ration p times then the superlattice will have a period 
D = di + d,2 + ■ ■ ■ + d n and L = Np layers. The inter- 
faces are perpendicular to the wave vector of the incident 
wave. The superlattice is placed between two media with 
the acoustic impedance Z medium and Z suhstrate . Meth- 
ods of calculation for the transmission and reflection co- 
efficients in a linear stratified medium are described in 
Refs. [§ H H P 
calculation in Fig 



We have presented results of our 
One should note that all filters 
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From the analysis in the preceding subsection it follows 
that the design of appropriate phonon filters is of primary 
importance. Using the analogy between reflection of a 
plane electromagnetic wave at the interface between two 
optical media with different refractive indexes n\ and 
and the acoustic wave reflecting at the boundary of two 
elastic media with different acoustic impendance Z\ and 
Z2, Narayanamurti et al. |lj, [ll| have developed appro- 
priate phonon filters for phonon spectrometery. 



FIG. 6: Band pass properties of designed filters (for the pa- 
rameters see Tab.fllJ). 



have an oscillatory behavior for the transmission coeffi- 
cient in the band pass region. This is typical for band 
pass filters. In Tab.|| we list impedance and thickness of 
the superlattices used in our calculations. Fig.|3| presents 
curves for the transmission spectra of the filter between 
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reservoir 2 and superconductor S with different acoustic 
impedance, layer thickness and layer number. Compar- 
ing with Fig.Q| we observe large phonon fluxes at fre- 
quencies about LObfq — A/h and uJbfq = 3A/H. The filter 
shown in Fig. [^d] may work well for those frequencies. 
The filter between reservoir 1 and superconductor S' is 
presented in Ref ]l7| |. The phonon filters for the SIN 
tunnel junctions can be designed by a similar analysis. 

TABLE I: Parameters of the presented niters (see Figure |^ 
). The impedance (Z) and the size of the layers (d) are given 
in units of 10 5 <?/(s cm 2 ) and in Angstroms (A), respectively. 





z x 


Z 2 


Zz 


z± 


^5 


di 


d 2 


d 3 


d4 


d 5 


N 


P 


a 


2.5 


1.8 


1.4 


0.9 


1.8 


50 


50 


50 


50 


30 


5 


11 


b 


3 


1.8 


1.3 


0.9 


1 


50 


50 


50 


50 


30 


5 


11 


c 


3 


1.8 


1.3 


0.9 


1 


40 


45 


50 


50 


30 


5 


11 


d 


3 


1.8 


1.3 


0.9 


1 


40 


40 


40 


45 


30 


5 


11 



B. Phonon filters and refrigeration 

The proposed PDE refrigerator is suited to refriger- 
ate an object attached to it. The effect of phonon filters 
can be estimated by the value of P/Pg from an experi- 
ment carried out by Pekola et al ||16|] where a membrane 
attached to the SIN tunnel junction was cooled. Our 
estimate for the net phonon absorption P/Pq by Equa- 
tion ([)]) varies from to 1CP 2 depending on temperatures 
T e and T p h- Table |l| presents our calculated values for 
the nonequilibrium tunnel junctions (see also Figure^). 
We find qualitative agreement between experiment and 
calculation for the value of absorbed power. 



TABLE III: Absorbed power P/P in the SIN tunnel junction 
at two different temperatures T — 0.5T C < and T = 0.8T C / , 
without and with phonon filter. The data below correspond 
to Fig. II . The ratio u/j= 0.02. 



T/T c , 


0.5 


0.8 


V/T c> 


1.31 


1.31 


effect 


total 


filtered 


total 


filtered 


P/Po 


-2.24 1Q-' Z 


-2.43 10" 2 


1.24 10" 4 


-2.5 W~ 2 



filtered phonon absorption for two different low pass fil- 
ters with transmittance 1 at ui < 2.95T C and u> < 3.25T C . 
Thus the phonon absorption is incresed in a wide tem- 
perature region. 




-0,05-1 



TABLE II: Maximal absorbed power P/P in the SIN tunnel 
junction at three different temperatures T — 0.27T C /,T = 
0.5T C / and T = 0.8T C / . 



T/T c , 


0.27 


0.5 


0.8 


V/T c , 


1.56 


1.31 


0.79 


v/l 


0.005 


0.02 


0.005 


0.02 


0.005 


0.02 


P/Po 


0.004 


0.014 


0.0052 


0.0224 


0.0011 


0.0042 



The phonon filters will increase the absorbed power 
by cutting off the emission window (see Figure ||) . Ta- 
ble III presents the absorbed power after placing a low 
pass filter (see Figures [| and between the SIN tunnel 
junction and an object. Phonon absorption is increased 
at T/T c = 0.5. At T/T c = 0.8 we see that SIN tunnel 
junction with phonon filter will absorb phonons (the mi- 
nus sign in Table III means absorption). Fig. m shows 
the temperature dependence of the total phonon absorp- 
tion P/Po for three different casses. The solid line is the 
phonon absorption by a SIN tunnel junction without 
phonon filter. The dotted and dashed curves present the 



FIG. 7: Phonon flux with low pass filter at lu 
uj = 3.25T C 



2.95T C and 



We recall that the refrigeration schema presented 
above is suited to refrigerate an object attached to it. But 
as follows from the phonon absorption spectra an analo- 
gous design with high pass filters can prevent the phonon 
absorption by SIN tunnel junction and thus reduce the 
heating of electrons by phonons (see Figures || 0). 



V. CONCLUSIONS 

The SSMR performance depends on various external 
and internal parameters. The PDE-based approach pre- 
sented here and the usual apprach [Q, [l6|, ^3) give the 
similar expression (Eq. ^|) for the net cooling in case of 
SIN tunnel junction. Thus one can conclude that the 
cooling of a membrane in the experiment by Pekola et 
al. pi is a direct implication of the phonon deficit effect. 
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The PDE— approach offers explicit spectral description 
for participating phonon fluxes, ft provides receipes to 
overcome unwanted heating channels and shows new de- 
tails which have been omitted in a simplified " cooling by 
extraction" description. 

We have analyzed the "phonon deficit effect" in SIN 
tunnel junctions. In contrast to SIS tunnel junctions 
there is a temperature region where the SIN tunnel 
junctions essentially absorb low frequency phonons. As 
temperature is increased, the SIN tunnel junctions also 
emit phonons to the environment. Depending on the 
value of the ratio u/'y there is net phonon emission or 
absorption. The temperature dependence of the net 
phonon absorption is consistent with experimental ob- 
servations [|| l||, there is heating at high and low 
temperatures and refrigeration at intermediate tempera- 
tures. 

To enhance the cooling power of the SIS(N) and tun- 
nel junctions and overcome the phonon heating at high 
temperatures, phonon filters can be effective. We per- 
formed calculations for the microrefrigerator phonon fil- 
ters. Those filters have superlattice design and can be 
obtained using contemporary nanotechnology methods. 
They can enhance the cooling power of SIS and SIN 
tunnel junctions. 



(compare with the formula (]?]), h = l,ks = 1) 
P = r u q po(0 P h - el (N° )du; q . 



(A.l) 



Here I ph el (N® ) is the phonon-electron coillision inte- 
gral and it is given by 



' Jo Jo 

(A.2) 



8e F 

f- 2S 2 5 (e - e' - u q )} , 



where Si and 52 are elementary collison processes. When 
the distribution function of electrons and holes is Fermi 
function then n e = n_ c . In this case Si and 52 can be 
written as 

5i(e, e', u>q) = 4 [(JV„ + l) n e n e , - N Uq (1 - n e ) (1 - n e .j\ 

(A.3) 



5 2 (e, e', oj q ) = 4 [(JV U + l) n e (1 - n c ,) - N„ q (1 - n e ) n e , 

(A.l! 



After substitutions of th e ex pressions (A.3), (A. 4) and 
(A.2) in the expression (A.l), the integration by delta 
function gives 
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APPENDIX: PHONON FLUX BETWEEN 
ELECTRONS AND PHONONS 



47tAwd 
8ef jo 



AiOqUlq Po(Wq) 



2S 2 (e,e- uj q ,LO q ) de 



Si(e,uj q - e,ojg) de 



(A.5) 



Finally evaluating these integrals by noting that T <C 
ud and using the integral representation of Rieman zeta 
function £(n + l)r(n + 1) = J °° x n / (exp(x) — l)dx one 
obtains 



In this appendix, we give the derivation of the expres- 
sion (||). 

The general expression for the phonon flux is given by 



P = 4 



A^r( 5 )c(5) (n-r; h ) 

lo7r epK 



(A.6) 
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